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Abstract 

The result provided in this paper helps complete a unified picture of the scaling 
behavior in heavy-tailed stochastic models for transmission of packet traffic 
on high-speed communication links. Popular models include infinite source 
Poisson models, models based on aggregated renewal sequences, and models 
built from aggregated on-off sources. The versions of these models with finite 
variance transmission rate share the following pattern: if the sources connect 
at a fast rate over time the cumulative statistical fluctuations are fractional 
Brownian motion, if the connection rate is slow the traffic fluctuations are 
described by a stable Levy process, while the limiting fluctuations for the 
intermediate scaling regime are given by fractional Poisson motion. 
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1 Introduction 



It is well-known that packet traffic on high-speed links exhibit data characteristics consistent 
with long-range dependence and self-similarity. To explain the possible mechanisms behind 
this behavior, various network traffic models have been developed where these features arise 
as heavy-tailed phenomena; see Resnick (2007) |15j . A natural basis for modeling such sys- 
tems, applied early on during these developments, is the view of packet traffic composed of a 
large number of aggregated streams where each source alternates between an active on-state 
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transmitting data and an inactive off-state. The traffic streams generate on average a given 
mean-rate traffic, they have stationary increments and they are considered statistically in- 
dependent. In particular, the transmission channel is able to accommodate peak-rate traffic 
corresponding to all sources being in the on-state. To capture in this model the strong positive 
dependence manifest in empirical trace data measurements, it is assumed that the duration 
of on-periods and/or off-periods are subject to heavy-tailed probability distributions. It is 
then interesting to analyze the workload of total traffic over time and understand the random 
fluctuations around its cumulative average. Our continued interest in these questions comes 
from the finding that several scaling regimes exist with disparate asymptotic limits. 

The first result of the type we have in mind is Taqqu, Willinger, Sherman (1997) [16J, 
which introduces a double limit technique. In this sequential scheme, if the on-off model is 
averaged first over the level of aggregation and then over time the resulting limit process is 
fractional Brownian motion. As the fundamental example of a Gaussian self-similar process 
with long-range dependence, this limit preserves the inherent long-range dependence of the 
original workload fluctuations. On the other hand, averaging first over time and then over 
the number of traffic sources the limit process is a stable Levy process. This alternative 
scaling limit is again self-similar but lacks long memory since the increments are independent. 
Moreover, having infinite variance the limiting workload is itself heavy-tailed. In Mikosh, 
Resnick, Rootzen, Stegeman (2002), [3], the double limits are replaced by a single scheme 
where instead the number of sources grows at a rate which is relative to time. Two limit 
regimes of fast growth and slow growth are identified and two limit results corresponding to 
these are established, where again fractional Brownian motion and stable Levy motion appear 
as scaled limit processes of the centered on-off workload. The purpose of this paper is to 
show that an additional limit process, fractional Poisson motion, arises under an intermediate 
scheme which can be viewed as a balanced scaling between slow and fast growth. In this case 
the scale of time grows essentially as a power function of the number of traffic sources. As 
will be recalled, fractional Poisson motion does indeed provide a bridge between fractional 
Brownian motion and stable Levy motion. 

The intermediate limit scheme discussed here is indicated in Kaj (2002) [8j, and introduced 
in Gaigalas and Kaj (2003) [7]j where limit results are given for a different but related class 
of traffic models under three scaling regimes referred to as slow, intermediate and fast con- 
nection rate. The workload process is again the superposition of independent traffic streams 
with stationary increments but now each source generates packets according to a finite mean 
renewal counting process with heavy-tailed interrenewal cycle lengths. The link to the class 
of on-off models is that each pair of an on-period and a successive off-period forms a renewal 
cycle and the number of such on-off cycles generate a heavy-tailed renewal counting process. 
Moreover, if we associate with each renewal cycle a reward given by the length of its on-period 
and apply a suitable interpretation of partial rewards, then the corresponding renewal-reward 
process coincides with the on-off workload process. 

To explain briefly the limit result in [7] under intermediate connection rate, let (N l (t))i 
be i.i.d. copies of a stationary renewal counting process associated with a sequence of inter- 



2 



renewal times of finite mean u and a regularly varying tail function F(t) ~ L{t)t~^ , charac- 
terized by an index 7, 1 < 7 < 2, and a slowly varying function L. Let m — > 00 and a — > 00 in 
such a way that mL(a)/a 7_1 — > uc 7 ' 1 for some constant c > 0. Then the weak convergence 
holds, 



1 t 1 

T(iV>i)-i)^-- C 7 7 ( i /c 



where Y^(t) is an almost surely continuous, positively skewed, non-Gaussian and non-stable 
random process, which is defined by a particular representation of the characteristic function 
of its finite-dimensional distributions. Additional properties of the limit process are obtained 
in Kaj (2005) [9 J and Gaigalas (2006) [6], where it is shown with two different methods that 
can be represented as a stochastic integral with respect to a Poisson measure N(dx, du) 
onRx M. + with intensity measure n(dx,du) = r ydxu~' y ~ 1 du. Indeed, 

Y i{ t )=i I hx,x+u]{y)dy{N(dx,du) - jdxu'^du), t>0, 
JrxR+ Jo 

([7] uses F(t) ~ L(t) r y~ 1 t~' y , consequently n(dx,du) = dxu~' y ~ l du). We call this process 
fractional Poisson motion with Hurst index H = (3 — 7)/2 6 (1/2, 1). With 

2 2 / \ 

CT 7 = ( 7 _i)( 2 -7)(3-7) = 2^(1-^(2^-1) = ^' (1) 

we may put Y~(t) = a^Puif) and obtain the standard fractional Poisson motion Pjj. A 
calculation reveals 

Cov(P H (s),P H (t)) = ^(\s\ 2H + \t\ 2H -\t- s\ 2H ). 
For comparison, fractional Brownian motion of index H has the representation 
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Bii(t) = — / / h x<x+u ](y)dyM(dx,du), 
&h JrxR+ Jo 

where M(dx, du) is a Gaussian random measure on M x M + which is characterized by the 
control measure (3 - 2H) dxu- 2i ~ 2 ~ H ^du. The covariance functions of E>h and Ph coin- 
cide. The fast connection rate limit for the model of aggregated renewal processes applies 
if mL{a)/a 1 ~ 1 — > 00 and the slow connection rate limit if mL(a)/a' y ~ 1 — > 0. For suitable 
normalizing sequences, the limit processes under these assumptions are fractional Brownian 
motion with Hurst index H = (3 — 7)/2 in the case of fast growth and a stable Levy process 
with self-similarity index I/7 in the slow growth situation, see [7|. 

A number of other models have been suggested for the flow of traffic in communication 
networks. The superposition of independent renewal-reward processes applies more generally 
to sources which attain random transmission rates at random times, and not merely switch 
between on and off. For a model where the length of a transmission cycle as well as the 
transmission rate during the cycle are allowed to be heavy-tailed, Levy and Taqqu (2000) 
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Pipiras and Taqqu (2000) [13| . and Pipiras, Taqqu and Levy (2004) |14| . established 
results for slow and fast growth scaling analogous to those for the on-off model. In addition, 
they obtained as a fast growth scaling limit a stable, self-similar process with stationary 
but not independent increments, coined the telecom process. A further category of models 
for network traffic with long-range dependence over time starts from the assumption that 
long-lived traffic sessions arrive according to a Poisson process. The sessions carry workload 
which is transmitted either at fixed rate, at a random rate throughout the session, or at a 
randomly varying rate over the session length. Such models, called infinite source Poisson 
models, are widely accepted as realistic workload processes for internet traffic. Indeed, it is 
natural to assume that web flows on a non-congested backbone link are initiated according 
to a Poisson process while the duration of sessions and transmission rates are highly variable. 
The conditions under which slow, intermediate and fast scaling results exist and fractional 
Brownian motion, fractional Poisson motion, stable Levy processes and telecom processes 
arise in the asymptotic limits are known in great detail for variants of the infinite source 
Poisson model, see Kaj and Taqqu (2008) [10J . In [10J, is called the intermediate telecom 
process. Mikosh and Samorodnitsky (2007) [5] consider scaling limits for a general class of 
input processes, which includes as special cases the models already mentioned as well as other 
cumulative cluster- type processes. It is shown that fractional Brownian motion is a robust 
limit for a variety of models under fast growth conditions, whereas the slow growth behavior 
is more variable with a number of different stable processes arising in the limit. 

Our current result completes the picture for the intermediate scaling regime, where neither 
of the mechanisms of fast or slow growth are predominant. In this case, where the system 
workload is under simultaneous influence of Gaussian and stable domains of attraction, we 
show that the fluctuations which build up in the on-off model are robust and again described 
by the fractional Poisson motion, parallel to what is known to be valid for infinite source 
Poisson and renewal-based traffic models. In the next section 2 we introduce properly both 
the on-off model and the renewal-based model to be used as an approximation and we state 
the relevant background results for these models. In section 3 we state the main result and 
give the structure of the proof. Section 4 is devoted to remaining and technical aspects of the 
proof. 



2 The on-off model and background results 

We begin by introducing the on-off model using similar notations as in [4] . Let X on , X\ , Xi , ■ ■ ■ 
be i.i.d. non-negative random variables with distribution F OQ representing the lengths of on- 
periods. Similarly let Y g, Yi,Y2, ■ ■ ■ be i.i.d. non-negative random variables with distribution 
F g representing the lengths of off-periods. The X- and Y-sequences are supposed to be 
independent. For any distribution function F we write F = 1 — F for the right tail. We fix 
two parameters, a on and a Q fr, such that 

1< a on < a off < 2, (2) 
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and assume that 

F on (x) = x- a ° n L on (x) and F oS {x) = x~ a ° s L oS (x), x -> oo, (3) 

with L on ,L Q ff arbitrary functions slowly varying at infinity. Hence both distributions F on and 
F ff have finite mean values /j, on and /i ff but their variances are infinite. Assumption ([2]) 
agrees with that of [4J. However, thanks to a simple symmetry argument, we can also cover 
the case a on > a g. The case a on = a ff , for which the on-off process is an alternating renewal 
process, falls outside of the class of processes we are able to study within the methodology 
developed here. 

We consider the renewal sequence generated by alternating on- and off-periods. For the 
purpose of stationarity we introduce random variables (Xq,Yq) representing the initial on- 
and off-periods as follows: let B, X i, Y^3 be independent random variables, independent of 
{X on , (X n ),Y g, (Y n )}, and such that B is Bernoulli with 

P(B = 1) = 1 - P(B = 0) = Mon //x, 

and Xon and have distribution functions 

F^(x) = / F on (s)ds and F^(x) = / F oS (s)ds, 

Mon 70 Moff JO 

respectively. Now, let 

X = BX£ and Y = BY oS + (1 - B)Y%. 

Note that Xq and Yq are conditionally independent given B but not independent. At time 
t = the system starts in the on-state if B = 1 and in the off-state if B = 0. With this 
initial distribution, the alternating renewal sequence is stationary and the probability that 
the system is in the on-state at any time t is /j, on / \i. Renewal events occur at the start of each 
on-period. Inter-renewal times are given by the independent sequence = Xi + Yi,i > 0, 
where Z{ has distribution F = F on * F Q g and mean fi = fj, on + // g for i > 1, and Zq has 
distribution function 

F e *(x) = - I" F(s)ds. 
V> Jo 

The renewal sequence (T n ) n >i with delay To is defined by 

n— 1 

T n = 'y Zi , 

i=0 

and we denote by N(t) the associated counting process 

N(t) = Y,Mo,t](Tn). 

n>0 



5 



Note that N(t) has stationary increments and expectation E[iV(i)] = t/fj,. Moreover, because 
of ([2]) , the tail behavior of the inter-renewal times is given by 

F(x) ~ L on (x)x~ a ° n , x ->• oo, (4) 

see Asmussen [I], Chapter IX, Corollary 1.11. The on-off input process is the indicator process 
for the on-state defined by 

I(t) = l [Q , Xo )(t) + £ l [T n ,T n+ X n+1 ){t), t > 0. 
ra>0 

The source is in the on-state if I(t) = 1 and in the off-state if I(t) = 0. The input process 
I(t) is strictly stationary with mean 

E[I(t)] = P(J(i) = 1) = fiojfi. 

The associated cumulative workload defined by 

W(t) = [ 1(a) ds, t>0 
Jo 

is a stationary increment process with mean E[Wt] = t/j, on / \x. 

Let (P , W 3 , N J )j>i denote i.i.d. copies of the input process /, the accumulative workload 
process W, and the renewal counting process N for the stationary on-off model. For m > 1, 
consider a server fed by m independent on-off sources. We define the cumulative workload of 
the m-server system as the superposition process 

rn 

W m (t) = Y,W j (t), t>0, m>l, 
i=i 

and the renewal-cycle counting process for m aggregated traffic sources by 

m 

N m (t) = ^N j {t), t>0, m>l. 
j=l 

In this paper, we are mainly concerned with the asymptotic properties of the cumulative 
workload when the number of sources, m, increases and time t is rescaled by a factor a > 0. 
Thus, we consider the centered and rescaled process 



W m (at) — mat /jL on / /j, 
b(a, m) 




t > 0, 



where the renormalization b(a, m) will be precised in the sequel. The asymptotic is considered 
when both m — > oo and a — > oo. The relative growth of m and a have a major impact on 
the limit. Let a = a m be the sequence governing the scaling of time and suppose a m — > oo as 
m — > oo (we will often omit the subscript m). Following the notation in [7J, we consider the 
following three scaling regimes: 
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• fast connection rate 

mL on (a)/a a ™- 1 -)■ oo; (FOR) 

• slow connection rate 

mL on (a)/a a ™- 1 ^ 0; (SCR) 

• intermediate connection rate 

mL on (a) /a a ° n ~ 1 -> ft c a ° n ~ 1 , < c < oo . (ICR) 

In [4], the asymptotic behavior of the cumulative total workload is investigated under 
conditions dFCRl) and (pCRl) . 

Theorem 2.1 (Mikosch a/.) Recall assumptions (dj) and |2|). 

• Under condition (jFCRp and wii/t i/te normalization b(a,m) = (a 3_a ° n L on (a)m) 1//2 , £/ie 
following weak convergence of processes holds in the space of continuous functions on 

W m (at) - matfi on /n /u on 
b{a,m) 

where Bn(t) is a standard fractional Brownian motion with index H = (3 — a on )/2. 

• Under condition (jSCRp and £/ie normalization 

b(a,m) = (l/F on )*"(am) := inf{x > : -F on (x) < 1/am}, 
then in the sense of convergence of finite dimensional distributions, 

W m (at) - matfi on / fi fdd [i E v , . n 

where X am ix(t) is a standard a OQ -stable Levy motion totally skewed to the right, i.e. 
such that 

*cw,i(l) ~S aoa (l,l,0), 

and 

aon _ r(2 - a on )cos(7ra on /2) 
1 - a on 

The intermediate regime for renewal processes was investigated in [TJ. The formulation 
adopted here is given in [9], and is an immediate consequence of (|4]). 

Theorem 2.2 (Gaigalas-Kaj) Under condition (jlCRp and with the normalization b(a,m) = 
a, the following convergence of processes holds: 

N m (at) - mat /n 1 

=> <r aon cP H {t c), 

a fi 



7 



where o" Qon is given in |ip and -Pff(i) is the standard fractional Poisson motion 

P H (t) = — [ f l [XtX+u] (y)dy(N(dx,du) -dxaonu-^^du) (5) 
°~a on JrxR+ Jo 

with Hurst index H = (3 — a on )/2. 



3 Intermediate limit for the on-off model 

In this section, we investigate the intermediate scaling limit for the on-off model. The following 
is our main result. 

Theorem 3.1 Under condition (jlCRp and with the normalization b(a,m) = a, the following 
convergence of processes holds in the space of continuous functions on M + : 

W m (at) - matij, on / n fj, oS 

o- aon cP H {t c), 

a fi 

with cr aon in |7P and Pn(t) the standard fractional Poisson motion in ([3p. 

Remarks 3.2 The fractional Poisson motion is not self-similar but does have a property 
of aggregate-similarity, introduced in [9], which allows for an interpretation of the scaling 
parameter c. Consider for each integer m > 1 the sequence c m = m. 1 ^" 011-1 ). Then 

m 

c m P H {t/c m ) f = Y, p h(t), 

i=l 

where Pjj,Pjj, . . . are i.i.d. copies of Ph. Consider also the sequence c' m = m~ 1 /( a ° n ~ 1 \ For 
any m, 

m 

J2^'mPh(t/c' m ) & P H {t). 

i=l 

Hence, by tracing the limit process in Theorem 13.11 as c m — > oo, we recover in distribution 
the succession of all aggregates Y2i<i<mPff> m — Also, by letting c' m —> we find that 
the limit process represents successively smaller fractions which sum up to recover fractional 
Poisson motion. 

These relations explain the fact that fractional Poisson motion acts as a bridge between the 
stable Levy process and fractional Brownian motion. First, {c H Pu(t/c)} converges weakly 

to {Bn{t)}, as c — > oo. Indeed, c^Pnit j c m ) ^= ^7= Yli<i< m P H^) anc ^ tne cen t ra l limit 
Theorem yields the Gaussian limit as m — > 00. The required tightness property is shown in 
[6|. Moreover, it is shown in [6] that c 1//a ° n P// (t/c) converges in distribution as c — > to the 
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o? on -stable Levy process. To see that the limit must be ct on - stable, take d = c ■ d m for any 
c> 0. Then 

c^Pnit/c) f = —^J2 dl/aonp H(t/d), m>l, 

1 = 1 

and, assuming that the rescaled process {c l l a ° u P}{{t/c))t>o converge to some non-trivial limit 
process L, we must have as c — > (and hence d — > 0) 

L{t) f = -^Vl'ft), m> 1. 

This indicates that the limit L must be a on -stable. 



Heuristics of the proof of Theorem 13.11 To motivate that the limit process in the inter- 
mediate connection rate limit appears naturally, we discuss a decomposition of the centered 
on-off process based on its representation as a renewal-reward model. We first note that the 
single source cumulative workload has the form 

N(t) 

W(t) =X /\t + J2 X i- ( T N(t)-i + X N{t) - t)+. 
i=i 

Similarly, focusing on off-periods rather than on-periods, we have 

N(t) 

t - W(t) = Y A t + Y, Yi - (T N{t) - t) A Y N(t) . 

i=l 

The centered single source workload is therefore 

W (t)-^t = -{t-W{t)) + ^t 
fj, /i 

N(t) 

= - / i off (iv(t)-t//x)-^(y i -/i off ) + i?(t) 

with 

R(t) = (T N[t) - t) A Y N{t) —Yq At. 
Thus, for the workload of m sources, 

m Ni(t) m 

W m (t) - f^mt = -Moff (iV m (i) - mtlii) EE W ~ ^ + E Rj ® ^ 

V j=l i=l j=l 

using obvious notations. The balancing of terms under the scaling relation (JlCR]), makes it 
plausible that both terms 

-E E 

j=i t=i j=i 
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vanish in the scaling limit. This suggests asymptotically, 

W m {at) - n on mat/n N m (at) - mat/fi 

~-Moff , ("J 

a a 

and so Theorem 12.21 would imply Theorem 13.11 In the next final section, we will compare 
rigorously the two processes in ([7]). 



4 Proof of Theorem 13.11 

The proof of Theorem 13.11 relies on the following three lemmas: 
Lemma 4.1 In the scaling (|ICR| . for all t > 0, 

EE c*7 o. (8) 

a j=l 1=1 

Lemma 4.2 In the scaling (lICRj) . for all t > 0, 

-YR j (at) =^0. 

Lemma 4.3 In i/ie scaling (jlCRp . i/ie sequence of processes 

W m (at) - mat/j, ori / n 



-, t > 0, m > 1 



a 

is in i/ie space of continuous functions on 



Proof of Theorem 13.11 By Theorem 12.21 Lemma [4. II and Lemma [4.2[ the convergence of 
finite-dimensional distributions , 

W m (at) - mat (j, on / n fdd Moff n / \ 

=> o- aon cP H {t/c), 

a n 

is a consequence of the decomposition given in (j6j) . By Lemma 14.31 the sequence is tight 
in the space of continuous functions on R + . Hence weak convergence holds in the space of 
continuous functions and Theorem 13.11 is proved. □ 
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Proof of Lemma 14.11 We construct an alternative representation of the random variable 
in the left hand side of ([8]) . Define 



AT 1 (at) = inf{fc > 0; + ^ Zj > at} 

i=l 

and for j > 2 

k 

N\at) = inf{A; > 0; x{ + £ Z± ^ > at}. 

i=l 

For m > 1, let N m (at) = YlJLi N 3 (at). The random variables N 3 (at),j > 1 are i.i.d and for 
each fixed t > 

m AP(at) iVm(ai) 
"EE ( y /-^ff) and - (^/-Moff), 

j'=l i=l j=l 

have the same distribution (note that the uni-dimensional marginal distributions are equal 
but not the multidimensional distributions). This representation will enable us to prove that 
under assumption (|ICR| . in the space of cad-lag functions on M + endowed with the Skorokhod 
topology, we have the convergence 

(amu \ 
£ - ) =>o. (9) 
i=1 ' u>0 

Moreover, 

— N m (at) => -. (10) 
am fj, 



Equations ([9]) and (|10p together imply 



iV m (at) 

- 53 (if - Moff ) 



a 



and this proves the lemma. Thus, it remains to prove (j9|) and (llOh . 

To this aim, recall that the random variables Y^,i > 1, are i.i.d. with distribution such 
that the tail function F g is regularly varying with index — a ff ■ Hence there exists a regularly 
varying function L such that the centered and rescaled sum 

amu 

X>/-Mofr) 



(am) l / a °«L{am) J u>o 



converges in the space of cad-lag functions to some a G fj-stable Levy process (see [12J, the exact 
form of L or of the limit process are not needed here). This implies the convergence property 
([9]), since a » (am) 1//Q ° ff L(am) under the scaling assumption (lICRj) with a on < aioff- 
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We now prove equation (jlOp . The stationary renewal process N(t) has mean t/fj, and 
variance given asymptotically by 



Var(JV(t)) ~ a 2 aon ^t 3 - a ° n L on {t), t -> oo 



see [TJ, Equation (30), and references therein. Hence, -^-N m {a£) has mean t//j, and variance 
under scaling (jlCRj) . such that 



Var 



N m (at) 

am 



Var [JV(at)] 



0. 



This shows that -^-N m {at) converges in distribution to t//j,, which is (jlOp . This ends the 
proof of Lemma 14.11 □ 



Proof of Lemma [4.21 Since - + (^ivim ~ *) A^jvj(t)> it is enough to prove that 

- m - m 

-E^o^O and -E( T lr|-*) Ay ^)^°- 

3=1 3=1 

By stationarity, the random variables 3^ and (T^- —t) A Y^,^ have the same distribution; 

they represent the remaining time after and t, respectively, of the first off-period after 0, 
and after t. Since both sums have the same distribution, we only consider the first one. 

Using Karamata's Theorem (see [3]), the tail function satisfies 
F^(x) = / F oS (s)ds -L oS (x) 

Moff Jx Moflf "off — 1 

as x —> oo. This implies that the random variable Yq has a regularly varying tail with index 
~{ a oS ~ 1) an d hence belongs to the domain of attraction of an (a g — l)-stable distribution. 
Therefore there exists a slowly varying function L, such that 

1 V Y j 

m i/Kff-i)L(m) frf 

converges in distribution to a stable law of index a Q s — 1 . Under scaling fllCRp with a on < a Q s , 
we have a » m l ^ a ° B ~ 1 ^ L{m) and hence 



7, T. Y « 



0. 
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Proof of Lemma 14.31 The proof given in [4] for fast scaling (IFCRP can be adapted to 
our settings. We recall only the main lines. According to Billingsley [2], Theorem 12.3, it 
is enough to prove that for any t\, t 2 with \t\ — t 2 \ < 1 and for some e > 0, there exists a 
constant C > and an ao > 0, such that for all a > ao 



E 



1 2 

- \{W m (at 2 ) - mai 2 /W V) ~ (Wm(a*l) - mati/j, on / fj,)\ 
a 



< c\t 2 -h 



l+e 



^Var[JU at ] < Ct 1+£ (11) 



Using the definition of W m , centering and stationarity, it is enough to prove that for all 
t E [0, 1] and a > ao, 

m, 

(the constant C may change from one appearance to another). However, according to [3], 
Equation (7.1), 

Vax(W t )~al m ?ft 3 - a °»L OJi (t), t -> oo. (12) 
This relation and the scaling (ICR) together imply, as a — > oo, 

~2 ~1 Cfcon ^1 C^on i .2 

_ „ £ a 3 " a -L on (a) ~ ^^-V Var[W a ], 

and so there is C > 0, such that for a large enough 

-^Var < C 

By (|12p . the function a i-> Var[W a ] is regularly varying with index 3 — a on - Then, using Potter 
bounds (see [3]), we conclude that there exist ao > and e < 1 — a on /2, such that for all 
t E (0, 1) and a > ao/t, 

V&r\W at ] 1 "in 
Var[W a ] ~ 1 - e 

(see the proof of Lemma 13 in [¥j for details). This implies that for all t E (0, 1) and all a 
such that at > oo, 

^Var[iy at ] < -SI— < 
a 2 1 — e 

On the other hand, if t < ao/a, then, for a large enough, 

™ VarfW ,1 < ° aH2 < C ^ < C 



and so 



m , „ a\- £ t 1+£ 



2 Var[H/ ai ] < C 2 _ q ° — < C t 1+£ . 
cr a z a ° n £ L on (a) 



In the last inequality, we use the fact that 2 — a on — e > and so a 2 ~ a ° n_£ L on (a) oo as 
a — 7- oo; taking ao large enough, we can suppose that for a > ao, a 2 ~ a ° n ~ £ L on (a) remains 
bounded away from zero. 

By combining the estimates for the cases at > ao and at < ao we obtain (fTT|) . which 
completes the proof. □ 
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